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Abstract 

The anomalous dimension matrix of dimensionally regularized four- 
quark operators is known to be affected by evanescent operators, which 
vanish in D = 4 dimensions. Their definition, however, is not unique, as 
one can always redefine them by adding a term proportional to {D — 4) 
times a physical operator. In the present paper we compare different defini- 
tions used in the literature and find that they correspond to different renor- 
malization schemes in the physical operator basis. The scheme transfor- 
mation formulae for the Wilson coefficients and the anomalous dimension 
matrix are derived in the next-to-leading order. We further investigate the 
proper treatment of evanescent operators in processes appearing at second 
order in the effective four-fermion interaction such as particle-antiparticle 
mixing, rare hadron decays or inclusive decays. 



1 Introduction 

In the past two decades much effort has been made to calculate QCD corrections 
to weak processes. The indispensable renormalization group improvement of 
perturbatively calculated Feynman amplitudes requires their factorization into 
Wilson coefficients and matrix elements, which are obtained from an effective field 
theory containing four-fermion interactions. When calculating QCD radiative 
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corrections to these four-fermion operators using dimensional regularization (D = 
4 — 2e) one faces evanescent Dirac structures such as 

7^7,7^(1 - 75) ® 7*7Y(1 - 75) - (4 + ae) 7AI (l - 75) ® 7 M (1 - 7s), (1) 

which vanish in D = 4 dimensions, but appear with a factor of 1/e in counterterms 
to physical operators. By introducing the parameter a in fll]) we have displayed the 
arbitrariness in the definition of the evanescent operators: A priori one can add 
any multiple of e times any physical operator to a given evanescent operator. In 
the literature one indeed finds different definitions of the latter. The consequences 
of this arbitrariness for renormalization group improved Green's functions are one 
of the main subjects of this paper. 

The role of evanescent operators in perturbation theory has been investigated 
since the pioneering era of dimensional regularization |I], 0, [|. When perturbative 
results are to be improved by means of the operator product expansion and renor- 
malization group (RG) techniques to sum large logarithms, new subtleties arise: 
First the matrix elements of evanescent operators can affect the matching equa- 
tion determining the Wilson coefficients which multiply the effective four-fermion 
operators |5|]. Second the appearance of evanescent operators in counterterms to 
physical operators and vice versa leads to the mixing of physical and evanescent 
operators during the RG evolution |2], |6| . In |2|, |5| a finite renormalization of the 
evanescent operators has been proposed to render their matrix elements zero. By 
this the Wilson coefficients of the evanescent operators become irrelevant at the 
matching scale. For this to be true at any scale it is important that simultane- 
ously the evanescent operators do not mix into physical ones. In || it has been 
proven for a very special definition of the bare evanescent operators that this 
is indeed the case, if the finite renormalization proposed in || [j| is performed. 
But does this feature hold for any definition of the evanescent operators, i.e. for 
any choice of a in ([!])? We will affirm this question in section |3| after setting up 
our notations and describing and generalizing the commonly used definitions of 
evanescent operators in section |[ 

It is well-known that a change in the renormalization prescription of the com- 
posite operators affects the Wilson coefficients and the anomalous dimension ma- 
trix in the next-to-leading order (NLO) and beyond. In section |] we will find that 
a change in the definition of the evanescent operators, i.e. a change of a in flip, 
leads to a different form of the physical part of the anomalous dimension matrix. 
Hence a different a corresponds to a different renormalization scheme. This result 
is of utmost practical importance for any calculation beyond leading logarithms, 
as it shows that it is meaningless to state a result for some anomalous dimension 
matrix without mentioning the definition of the evanescent operators used in the 
calculation. If one wants to combine some anomalous dimension matrix with 
Wilson coefficients or perturbative matrix elements calculated with a different 
definition of the evanescent operators, one clearly needs scheme transformation 
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formulae for the Wilson coefficients and the anomalous dimension matrix. We 
will derive these scheme transformation formulae in the next-to-leading order in 
section |, too. 

When studying particle-antiparticle mixing or rare decays one faces Green's 
functions with two insertions of four-fermion operators. The second main subject 
of this paper is to work out the correct treatment of these Green's functions when 
one or both inserted operators are evanescent. In section |^ we extend the results 
of || m and of sections ^| and |] to this case of double insertions. 

Then in section ^| inclusive decays are discussed. We close our paper with our 
conclusions. 



2 Preliminaries and Notation 

Let {Qk = ipqk^P • 4>qk4>, k = 1,2,3,...} be a set of physical dimension-six 
four-quark operators. We are interested in the Green's functions of a SU(N) 
gauge theory with insertions of Qk renormalized by minimal subtraction (MS). 
The arguments are easily generalized to other mass-independent renormalization 
schemes like MS. The Dirac structures Qk = % ® % corresponding to Qk are 
considered to form a basis of the space of Lorentz singlets and pseudosinglets 
for D = 4. Neither the Lorentz indices of qk and q~k are displayed nor any 
flavour or colour indices, which are irrelevant for the discussion of the subject, 
[r ® 1] Qk [1 ® r'] means Tqk <S> <?aT' '. Frequently we will use the example 

Q = T/i (1-75)® 7"(l-7 5 ). (2) 

The matrix elements of Qk have some perturbative expansion in the gauge 
coupling g: 

zi(Q h r) = E(if^y<or> w - ( 3 ) 

Here is the quark wave function renormalization constant. The right hand side 
of (|3]) still contains divergences, which are to be removed by the renormalization 
of the operators Qk ||. 
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Now the insertion of Qk into the one-loop diagrams of fig. [1] yields a lin- 
ear combination of the Qi's and a new operator with the Dirac structure 

Q' k = [i P ia®i}Q k [l®rY}: 

(Q h r) (1) = 4?(4) (0) +<^(^) (0) no sum on k, (4) 

where (. . denote tree level matrix elements. Both coefficients have a term 
proportional to 1/e and a finite part. Q' k is now decomposed into a linear com- 
bination of the Qi's and an evanescent operator: 

<?r c = (fki + ak^Q^ + EilQ^ + Ois 2 ). (5) 

Here the /^'s are uniquely determined by the Dirac basis decomposition in D = 4 
dimensions. The a^'s, however, are arbitrary, and a different choice for the djy's 
corresponds to a different definition of Ei[Q k ] = E%[Qk, {a r s}]- When going be- 
yond the one-loop order new evanescent operators E 2 [Qk], E 3 [Qk], . . . will appear. 
Their precise definition is irrelevant for the moment and will be given after fll6|) . 

Now in the framework of dimensional regularization the renormalization of 
some physical operator Q k requires counterterms proportional to physical and 
evanescent operators: We define the renormalization matrix Z by 

ZUQ h D = Z M {Qr) + Z^jEjlQ^™) . (6) 

Here and in the following we will distinguish the renormalization constants related 
to some evanescent operator Ej[Q m ] by denoting the corresponding index with 
Ejm- (ED and @ imply that depends on the a rs 's, while is independent 

of them. We define the coefficients in the expansion of Z in terms of the gauge 
coupling constant g and in terms of e by 




The first analysis of evanescent operators in the context of RG improved 
QCD corrections to electroweak transitions has been done by Buras and Weisz 
[f|. They have determined the a^s by choosing some set of Dirac structures 
M = {'j^ . . . ®7 ( - 10 ' ) }, which forms a basis for D = 4, and contracting 

all elements in M with Q' k and Qi in (||): 

tr ( 7 (m) ^7 M <?l) = (/« + a kl e) tr ( 7 (ni) «7 (m) «) + 0(e 2 ), 

no sum on k and on m—1,. . . , 10. (8) 

The solution of the equations (|8|) uniquely defines the j^i + a^e. In other words, 
Ei[Qk) obeys the equations: 

ExiQkhrslif^ = 0(e 2 ) for m=l, ...,10, (9) 
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where i,j, r, s are Dirac indices. 

Our arguments will not depend on the scheme used for the treatment of 75. 
In the examples we will use a totally anticommuting 75. This does not cause any 
ambiguity in the trace operation in @, because all Lorentz indices are contracted, 
so that the traced Dirac string is a linear combination of 75 and the unit matrix. 

E.g. the choice of 

M = {1® 1,1® 75, 7s® 1,75®75,7a»®7 a *>7m®7 m 75, 



7s7m » T, 757^ » 7"7s, <V ® ^ 1 75<V ® ( 10 ) 
gives for Q in @ 

Q' = (1 - 75) ® 7 M 7V (1 - 75) = (4 - 8e)Q + ^[Q] + 0(e 2 ) (11) 

as in ||. We remark that this choice a = —8 respects the Fierz symmetry, which 
relates the first to the second diagram in fig. |l[ 

A basis different from M in (|T0|) yields the same fk/s, but different a^'s. For 
example by replacing the sixth and eighth element of M in fllCf ) by ^ a li3l5®l a l l3 l S 
and 757q7/?75 ® ^ a ^^ & one finds 

Q' = 4Q + 16e(l + 7 5)®(l-75) + ^[Q fc ] + 0(£ 2 ), 

instead of (|TTJ) , i.e. a different evanescent operator. The Dirac algebra is infi- 
nite dimensional for non-integer D and is spanned by M and an infinite set of 
evanescent Dirac structures. Hence one can reverse the above procedure and first 
arbitrarily choose the a^s and then add properly adjusted linear combinations 
of the evanescent structures to the elements of M such as to obtain the chosen 

Yet the so defined evanescent operators do not decouple from the physics in 
four dimensions: In |fj, |S| it has been observed that their one-loop matrix elements 
generally have nonvanishing components proportional to the physical operators 

Qk'- 



(£i[Q fc ] barc > 



bare \ (1) 



(^) (0) + - z\ 



1 



+ - 

e 



(i) 



(E 2 [Qk 



+ 0(e). 



(Ei[Qk\) 



(o) 



(12) 



Here a second evanescent operator E 2 , which will be discussed in a moment, has 
appeared. Clearly no sum on k is understood in (O) and in following analogous 
places. In ( |T2"| ) [Zq ]E lk ,i is local, because it originates from the local 1/e-pole of 
the tensor integrals and a term proportional to e stemming from the evanescent 
Dirac algebra. For the same reason there is no divergence in the term proportional 
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to {Qi) {0) . Now in @,§ it has been proposed to renormalize E\ by a finite amount 
to cancel this component: 



E 1 [Q k ] h -+ ' ! 



16tt- 



1 r 

£ 
1 
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(1) 



(13) 



With ( |i~3"|) the renormalized matrix elements of the evanescent operators are 0(e), 
so that they do not contribute to the one-loop matching of some Green's function 
G rcn in the full renormalizable theory with matrix elements in the effective theory: 



iG™ = C l (Q l )- + C Elk (E 1 [Q k ])- + 0(g 4 ) 



(14) 



i.e. the coefficients C# lfc are irrelevant, because they multiply matrix elements 
which vanish for D = 4. In [||] it has been further noticed that in (pf) 
influences the two-loop anomalous dimension matrix of the physical operators, 
so that the presence of evanescent operators indeed has an impact on physical 
observables. In a different context this has also been observed in [[7). 

Next we discuss E 2 [Qk], which has entered the scene in (|12|): When inserting 
Ei[Qk) defined in (j5|) into the one-loop diagrams of fig. [I], one involves 



Ok 



[7p7.®1] Q'k [1®tV] 

[/ + ae\ 2 H Qi + (f H + a H e) E x [Qi] 

+ [7 p7(T ®l] E l [Q k ] [l® 7 CT 7l 

{[/ + + hie} Q t + E 2 [Qk] + O (V 



(15) 
(16) 



which defines E 2 [Qk] = E 2 [Qk, {a rs }, {Ks}]- Only the last term in fll5[) can 
contribute to the new coefficients bki- If the projection is performed with e.g. M 
defined in fllCf) , one finds bki = 0. In our discussion we will keep bki arbitrary. 
Clearly, one has a priori to deal with the mixing of an infinite set of evanescent 
operators {^[Qfc]} for each physical operator Qk, where Ej + i[Qk] denotes the 

new evanescent operator appearing first in the one-loop matrix elements of Ej [Qk] ■ 
With the finite renormalization of Ei[Qk] in ( |i~3"D the evanescent operators 
do not affect the physics at the matching scale, at which ( |T4] ) holds. In or- 
der that this will be true at any scale /x, however, one must also ensure that 
the evanescent operators do not mix into the physical ones. This has been 
noticed first by Dugan and Grinstein in [0]. For the operator basis (Q,E) = 

y(rn) g m the quantity 



^This is the case for any basis M in which for each 



(m) 



i 7 



is a linear combination of the elements in M. 
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(Qi, . . . Q n , Ei [Qi] , . . . Ej [Qk] , • • •) this means that the anomalous dimension 
matrix 
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1QQ 1QE 
lEQ JEE 



has an upper block-triangular form with ^eq — 0. 

The authors of |J have introduced another way to define the evanescent op- 
erators, which is also frequently used: It is easy to see that one can restrict the 
operator basis {Qk} to the set of operators whose Dirac structures qk, qk are com- 
pletely antisymmetric in their Lorentz indices. This is the normal form of Dirac 
strings introduced in ]3[]. Dirac strings being antisymmetric in more than four 
indices vanish in four dimensions and are therefore evanescent. Operators with 
five antisymmetrized indices correspond to E\ in our notation, and E2 would be 
expressed in terms of a linear combination of Dirac structures with seven and 
with five antisymmetrized indices. Clearly this method also corresponds to some 
special choice for the a^'s and fe^'s in (|5|) and (|i6|). Now in || the authors have 
proven that with the use of those definitions and a finite renormalization analogous 
to ( p~3| ) the anomalous dimension matrix indeed has the desired block-triangular 
form, so that the evanescent operators do not mix into the physical ones. While 
the anomalous dimension matrix is trivially block-triangular at one-loop level, 
the proof for the two-loop level was given in || by the use of the abovementioned 
special definition of the evanescent operators. The latter, however, has some very 
special features, which are absent for the general case with arbitrary a^s and 
bki's, e.g. the definition used in || automatically yields an anomalous dimension 
matrix which is tridiagonal in the evanescent sector. 

Consider now a definition of the evanescent operators different from the one 
used in 0: By inserting the definition (^) of Ei [Qk] into (0) one realizes that 
depends on the a^'s. Similarly at the two-loop level Z\ depends 

E lk ,l L J Ejk,l 



^0 



on the definition (^), so that one has to wonder which choices for the a^'s lead 
to the desired block-triangular form of 7 with ^eq — 0. In the following section 
we will prove that any choice is permissible. Further we will find that also the 
bkis may be chosen completely arbitrary. 

On the other hand the physical submatrix 7qq of the anomalous dimension 
matrix depends on the a^'s as we will show in section |j. Hence the freedom in 
the definition of the bare evanescent operators induces a renormalization scheme 
dependence in the physical sector of the operator basis. This feature has not 
been discussed in the literature so far. As emphasized in the introduction it is 
of practical importance for NLO calculations to know the scheme transformation 
formulae for the physical submatrix jqq and the Wilson coefficients. We will 
come back to this point in section f|. 
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3 Block Triangular Anomalous Dimension Ma- 
trix 



Consider some set of physical operators {Qk} which closes under renormalization 
together with the corresponding evanescent operators {Ej[Qk\ '■ j > 1}- Their 
0(£)-parts chosen arbitrarily. We want to show that the block of 

the anomalous dimension matrix describing the mixing of Ej[Qk] into Qi equals 



zero, 



'jkr 



0. 



(17) 



provided one uses the finite renormalization described in (|13|). 

Our sketch will follow the outline of ||, where (|Tfl) has been proven by com- 
plete induction. At the one-loop level ( ff7|) is trivial, and the induction starts in 
two-loop order: The next-to-leading order contribution to the anomalous dimen- 



sion matrix 7 



16tt 



7 



(i) 



(0) 



+( 



I671- 2 



7 (D + 



reads 0: 



-4Zf)-2/3 4 1) + 2{z«,4 1) } 



The nonzero contributions to (|i~7|) in two-loop order are 
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(i) 



-4 



- J jki lj 



(2) 



Ej k ,l 



-2/? 



Ej k ,l 



+2 Z\ 



r(l) 



E rs ,l 



+ 



E ik ,m 



(1) 



ml 



(19) 



Here (p!9|) contains terms which are absent when the special definition of the 



is used: In 



one has 



for j > 2 



evanescent operators m 

contrary to the general case, where any evanescent operator can have counterterms 
proportional to physical operators. 



Next we look at 



(2) 



, which stems from the 1 / e-term of the O ( g 4 )-matrix 

Ej k ,l 

elements of Ej[Q k }. As discussed in M, these 1/e-terms originate from l/e 2 -poles 
in the tensor integrals multiplying a factor proportional to e stemming from the 
evanescent Dirac algebra. Now in each two-loop diagram the former are related 
to the corresponding one- loop counterterm diagrams by a factor of 1/2, because 
the non-local 1/e-poles cancel in their sum ||. For this to hold it is crucial 
that the one-loop counterterms are properly adjusted, i.e. that they cancel the 
1/e-poles in the one- loop tensor integrals. In the one- loop matrix elements of 
evanescent operators the latter are multiplied with e originating from the Dirac 
algebra. Hence the proper one-loop renormalization of the evanescent operators 
must be such as to give matrix elements of order s, as shown for Ei[Qk] in (JT5|). 
From the one-loop counterterm graphs one simply reads off: 



r(2) 



-(1) 



E jk ,m 



(1) 



ml 



z\ 



(1) 



Ej k ,E r 



AD 



E rs ,l 



-Pa 



Ej k ,l 
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which yields the desired result when inserted into ([19]). Here the first two terms 
stem from insertions of physical and evanescent counterterms to Ej[Qk], while the 
term involving the coefficient of the one-loop /3-function (3(g) = —g 3 /(16ir 2 )f3o 
originates from the diagrams with coupling constant counterterms. The terms 
involving the wave function renormalization constants cancel with those stemming 
from the factor Z\ in (|6]). 

The inductive step in |J proving (fL7D to any loop order does not use any 
special definition of the evanescent and therefore applies unchanged here. 



4 Evanescent Scheme Dependences 



In this section we will analyze the dependence of the physical part of j^ 1 ' given 
in (|i"8|) and of the one- loop Wilson coefficients on an and bu. In practical next- 
to-leading order calculations one often has to combine Wilson coefficients and 
anomalous dimension matrices obtained with different definitions of the evanes- 
cent operators and it is therefore important to have formulae allowing to switch 
between them (see e.g. appendix B of @). 

We start with the investigation of the dependence of on an- The bare 
one-loop matrix element 



(i) 



+0(e) 



(i) 



d, 



(i) 



(Qj) 



(0) 



(i) 



k.E, 



mQ k ]) {o) 



(20) 



is independent of an, which is evident from (Q). Ei[Qk] depends linearly on an 
through its definition (H) with the coefficient 



_d_ 

dai 



A[Q h 



-e 5ki Qi 



so that 



gives: 



_d_ 

da L 



4 



(i) 



Z\ 



(i) 



(21) 



(22) 



while Zi is independent of an. 

In the same way on can obtain the a^-dependence of Z\ . @ reads to 
two- loop order (cf. @): 

(Q h n {2) = 4?(4) (0) + 4lJ^i[Q m ]) (0) + 4lj^[g m ]) (0) 

+zt ) (Q7) (1) + Z^jE^rY-V + (QT) (2) - (23) 

From (|I^) we know 

d 



da, 



-E 2 [Q r 



~£ [fmi^lj &miflj\ Qj; 



(24) 
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and from ( |20| ) one reads off: 

(<T> (1) 



d, 



(i) 



(Qj) 



(o) 



(25) 



These relations and (p2|) allow to calculate the derivative of (^) with respect to 
an. Keeping in mind that the evanescent matrix elements are 0(e) the 0(l/e)- 
part of the derivative yields: 



d 
dan 



(2) 



kj 



+ 



J ki L 



zf 



zf 



k,E u 



Si ; 



k,E 2r 



(Smifij + fmAj) , no sum on i. (26) 



Again 



,(2) 



can be extracted from the one-loop counterterm diagrams as 



described in the preceding section: 



4 2) 



,(2) 



k,E, 



z[ 1] ] \z? 

J ki L 



1 r 



z\ 



(1) 



z\ 



(1) 



2#> Z\ 



(i) 



no sum on i 



Eu, En 



>ki 



k,E 2r , 



z\ 



(1) 



k,E, 



z\ 



(1) 



Eik,E 2 k 



5, 



km: 



no sum on k. 



After inserting (|27|) into (£6|) we want to substitute the last term in 
this we derive both sides of (|l^) with respect to an giving: 



(27) 
For 



Z\ 



(i) 



Elk,E 2k 

d 



(Skifij + fki^ij) 



da. 



ii 



(i) 



+ 



(1) 



Li 



6, 



hi 



z\ 



(1) 



SkiSij, no sum on k. (28) 



Finally one has to insert the expression for (|26|) obtained by the described sub- 
stitutions into 



_d_ 

dai 



7 



(i) 



-4- 



Z] 



(2) 



+ 2 



(i) 



g 

fc, E ife dan 



no sum on k. 



which follows from (|TH). The result reads: 



_d_ 

dai, 
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(i) 



ki 



-2 



Z\ 



(i) 



^(1) 



+ 2 z} 1 Z{ 



(1) 



i,Eu 



+2/3 Zi 



(i) 



i,E u 



no sum on i. 



(29) 



Since the quantities on the right hand side of (|29| ) do not depend on an, one can 
easily integrate (|9|) to find the desired relation between two 7's corresponding 
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to different choices for au in ([5]). To write the result in matrix form we recall the 
expression for the physical one-loop anomalous dimension matrix 



(0) 



7 



and introduce the diagonal matrix D with 



(i) 



D, 



(i) 



i,Ei 



no sum on %. 



Hence 



7 W(o' 



7 (1) (a)+ \D- (a'-a),7 (0) l + 2 (3 D ■ (a' - a), 



(30) 



(31) 



where the summation in the row and column indices only runs over the physical 
submatrices. 

( |31~D exhibits the familiar structure of the scheme dependence of [ 10(1 . Usu- 
ally scheme dependences are analyzed for a fixed definition of the bare operators 
and different subtraction procedures. Our situation, however, is more compli- 
cated, because we investigate the scheme dependence associated with different 
definitions of the bare operator basis (i.e. of the bare evanescent operators). 

The dependence of the one-loop matrix elements on a can be found easily 
from (E5D and (E3): 
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1 + ^—D ■ (a' - a) 



{Q)-(a) + 0{g*). 



(32) 



Since in ([TJP G does not depend on a and the evanescent matrix element is 0(e), 
the corresponding relation for the Wilson coefficients at the matching scale reads: 



C T (a') = C T (a) 



16vr 2 



:D ■ (a' 



+ 0(g 4 ). 



(33) 



Hence we can apply the result of [10], which shows that in the renormalization 
group improved Wilson coefficient the scheme dependences in ([JT]) cancels the 
one in ([32]), so that physical observables are scheme independent, provided the 
hadronic matrix elements are defined scheme independently. 

Let us drop some words on the results ([H]), ([32]) and (|33|): In general one 
would expect scheme transformation formulae involving the full operator basis 
(Q, E). Yet all summations only run over the indices corresponding to the physical 
operators, the only ingredient from the evanescent sector being the matrix D. 
This is why we could not simply deduce ([Jl]) from (|32| ) using the results of [|1C 



Possible contributions from summations over evanescent operator indices in the 
matrix products in ( |3~Tp cannot be inferred from (^), because there they would 
multiply vanishing matrix elements. 
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In the same way one can investigate the dependence of 7 < - 1 - ) on bu given in 
fllpp: While z[ 2 ^ and depend on bu, this dependence cancels in (^8j). Hence 



a- 



neither -y ( ' nor the one-loop Wilson coefficient are affected by the choice of b. 

In general and do not commute, so that one has to cope with com- 
plicated matrix equations in order to solve the renormalization group equation 
in next-to- leading order [|I0 |. Now one can use fl3T|) to simplify 7^: By going to 



the diagonal basis for 7^ = diag (7* ) one can easily find solutions for a! — a in 
( |31~D which even give 7 < - 1 - ) (a / ) = provided that all Z^e^s are nonzero and all 
eigenvalues of 7*- - 1 satisfy 7}°^ — 7 J -°' ) 7^ 2(3$. We will exemplify this in a moment. 

A choice for an which leads to a 7^ commuting with 7^°-' has been done 
implicitly in ||: There the mixing of the two operators Q + — Q (l + lj /2 and 

Q_ = Q (l - l) /2 has been considered, where 1 and 1 denote colour singlet and 
antisinglet and Q was introduced in @. Now Q + is self-conjugate under the Fierz 
transformation, while Q_ is anti-self-conjugate, so that 7^ is diagonal to maintain 
the Fierz symmetry in the leading order renormalization group evolution. As 
remarked after (|TT|), the definition of in (|TT|) is necessary to ensure the Fierz 

symmetry in the one-loop matrix elements. Consequently with (JTTJ) also 7^^ has 
to obey the Fierz symmetry preventing the mixing of Q + and Q_, i.e. yielding a 
diagonal 7^ 1 - ) . A different definition of £i[Q] would result in non-Fierz-symmetric 
matrix elements, but in renormalization scheme independent expressions they 
would combine with a non-diagonal 7M such as to restore Fierz symmetry. 

Let us consider the example above to demonstrate that one can pick a' such 
that 7 ( - 1 )(a') = 0: In || the definitions 

Ei[Q±] = (±l-^)[ 7p 7 CT 7,(l-75)®7V7 P (l-7 5 ) 
-(4-8e) 7,(l-75)®7"(l-75)], 

i.e. a ++ = 8(1/ N — 1), a = 8(1/ N + 1), a + _ = a_ + = 0, were adopted yielding 

a diagonal 'y^(a) =diag(Vj!' ) (a), ^^(a)^ . From the insertion of Q + and Q_ 
into the diagrams of fig. |l| one finds Z +tEl+ = Z_ El _ = 1/4. Hence if we pick 



a 



±± = a±± - 2/7^ (a) /p and a' ± = 0, (|T|) will imply 7 (1) (a') = 



5 Double Insertions 



5.1 Motivation 



In the following we will investigate Green functions with two insertions of local 
operators. Consider first the effective Lagrangian to first order written in terms 
of bare local operators 



C l - C k Z k iQ\" c + C k Z k E rl E r [Qi 



bare 
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+ c E ]k ZJ l Q h l " c + C Ejk Z E Erl E r [Qi 



(34) 



According to the procedure presented in the preceding sections, the coefficients 
C E ., were found to be irrelevant and therefore remained undetermined. 

Now consider 4-fermion Green functions with insertion of two local operators 
from C l 



T^ 2 {^Jd D yC l (x)C l (j/))* 3 *4 



. 



(35) 



Such Green functions appear in applications like particle-antiparticle mixing or 
rare hadron decays. The diagram contributing to lowest order is depicted in 
fig. 0. Renormalization of them in general requires additional counterterms pro- 




Figure 2: The lowest order diagram contributing to the Green function in fl35"l) . 

portional to new local dimension-eight operators Qi, because the diagram of fig. | 
is divergent: 



C = C l + C 11 



(36) 



Qi 



C n J ^—1 *Obarc 

+C k C E ., k , ^Z k ^ iy ffi™ c + Z k ^ iki Eri E r Q 
+C E]k C Ejiy \Zz l ]kE]iki l Q h ^ c + Z E * kE . fkhEri E r [Q 
+C k Z kl Q z arc + C k Z kEr E r yQi 
+ c E ]k Z E ] kl QY IC + C E]k Z E l 3kEr E r [Q 



(37) 



Here Z~ X Q. are the local operator counterterms needed to renormalize the diver- 
gences originating purely from the double insertion. Further we have explicitly 
distinguished physical and evanescent operators. The renormalization constants 
Z_, r , clearly being symmetric in their first two indices, give rise to an inhomo- 
geneity in the RG equation for the Wilson coefficients C k , C Ers , which we call the 
anomalous dimension tensor of the double insertion. Note, that this quantity also 
has three indices, see (f4"E|). It has become standard to define the local operator Qi 
with inverse powers of the coupling constant such that Z~ l = O (g 2 ) to avoid mix- 

ing already at the tree level. As an example take Q — 7 M (1 — 75) <8>7 M (1 — 75) 
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for which G\ = (g 2 ). For simplicity, we assume the QiS to be linearly indepen- 
dent from the Q fc 's 0- The E r Qi in ([37|) are defined analogously to fl5|) with new 

coefficients fki, au, bki, etc. Hence new arbitrary constants dki, bki potentially 
causing scheme dependences enter the scene. 
Clearly the following questions arise here: 

1. Are the coefficient functions G-E ih irrelevant also for the double insertions; 
i.e. do 

' TEQ\ and Q TEE^j (38) 

contribute to the matching procedure and the operator mixing? 

2. Does one need a finite renormalization in the evanescent sector of double 
insertions; if yes, how does this affect the anomalous dimension tensor? 

3. How do the C\ and anomalous dimension matrices depend on the aki, bki, 
aw, bki ? 

4. Are the RG improved observables scheme independent? 
5.2 Scheme Consistency 

In this section we will carry out the program of section |] for the case of double 
insertions to answer questions [l] and ||] (on page |14"|) . 



Two cases have to be distinguished: The matrix element of the double insertion 
of the two local renormalized operators can be divergent or finite: 

1 f Tn n \ I diver g ent > case 1 /o Q \ 

2 J TQkQk ') = { finite , case 2 " (39) 

Case 1 is the generic one, appearing in the calculation of the coefficient 773 in 
K°-K° mixing or in K — > -nvv [1^. Case 2 appears, if the divergent parts of 



different contributions to (|39D add such that the divergences cancel. It is realized 
e.g. in the determination of 771 in K°-K° mixing [III]. Therefore we need or do 



not need a separate renormalization for the double insertion 

, case 1 ^ 
, case 2 

Since we need an extra renormalization in case 1, let us introduce the symbol 
/ T QQ for the completely renormalized operator product constructed from 




^e.g. the Qfc' s represent AF = 1 operators, the Q;'s denote AF = 2 operators, where F is 
some quantum number, which is conserved by the SU(N) interaction. 
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two renormalized local operators Q with an additional renormalization factor for 
the double insertion. 

Let us start the discussion with the matching procedure: At some renormal- 
ization scale we have to match Green functions obtained in the full theory with 
Green functions calculated in the effective theory: 



iG 1 ' 



CkCk 

+c Ejl (e s \q 



T QkQw 



T Ei [Qk] E v [Q h 



T QkEi' [Qk' 
Q (Qi . 



(41) 



where G rcn corresponds e.g. to a "box" function in the full SM. Since the coeffi- 
cients CE jk must be irrelevant for this matching procedure, one must have 

O (e°) in case 1 (LO) 



T Ej [Qk] Qk 



0{e l ) 
0(e l \ 



in case 1 (NLO and higher) 
in case 2 



(42) 



and analogously for two insertions of evanescent operators. To understand this 
recall that the purpose of RG improved perturbation theory is to sum logarithms. 
In case 1 the LO matching is performed by the comparison of the coefficients of 
logarithms of the full theory amplitude and the effective theory matrix element 
([42]) (the latter being trivially related to the coefficient of the divergence), while 
the NLO matching is obtained from the finite part and also involves the matrix 
elements of the local operators flT2L [T3|. In case 2 the matching is performed with 
the finite parts in all orders [ITT] . In both cases the condition (|42"D is trivially ful- 



fil both cases the condition ([42]) 
filled in LO, because the evanescent Dirac algebra gives an additional e compared 
to the case of the insertion of two physical operators. Therefore a finite renormal- 
ization for the double insertion turns out to be unnecessary at the LO level. This 
statement remains valid at the NLO level only in case 2, in case 1 condition 
no longer holds if one only subtracts the divergent terms in the matrix elements 
containing a double insertion. With the argumentation preceding ( p"3[) one finds 
that in this case the finite term needed to satisfy the condition (|42"D is local and 
therefore can be provided by a finite counterterm. 

The operator mixing is more complicated. To deal with this, we need the 
evolution equation for the Wilson coefficient functions Ck, Ce ts , which can be 
easily derived from the renormalization group invariance of C u and reads 



d/x 



with the anomalous dimension tensor of the double insertion 



lkn,l = — [ikk'Snn 1 + $kk'lnn'] Z k , n , yZyi — 



d 7 



-1 
knl' 



Zv 



/'/• 



(43) 



(44) 
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Using the perturbative expansions for the renormalization constants 



^knl 



E 



-l,(j) 



16tt 2 



J knl 



J kn,, 



j 1 

E- 



-i,Cj)1 



and Z we derive the perturbative expression for 



knl 



(45) 



7fenJ 



167T 



2 7S, + 



16vr 2 



7fcnJ + • • • 



(46) 



in (M) up to NLO: 



/ kn,l 

J 1 ) 
lkn,l 



•1.(1) 



-1.(2) 



knl 



kn,l 
-1.(1) 



-2 [z^* 1 
-2{[^ (1) 

-2{[zr 1 - (1) 



+ 2 £ 



-1.(1) 



knl 



-1,(1) 



fefc' 
kk 



J 

£-1,(1) 
1 JZ'Z 



knl 



■1,(1) 



fcn,F 



-1,(1) 





1,(1)" 






1,(1)" 




tin' J 






1,(1)" 


nn' 1 


^0 


1,(1)" 



fc'n'i 



(47) 



The indices run over both physical and evanescent operators. The reader may 
have noticed, that we have used the perturbative expansions of Z* -1 , Z^ 1 rather 
than Z, Z as in the previous sections. This is more convenient for the case 
of double insertions. Using these equations, the finite renormalization ensuring 
(f42|) to hold and the locality of counterterms, one shows in complete analogy to 
section H: 



" E r fJ,,n 







and 



I E-uE.i ,n 



0. 



(4* 



i.e. a double insertion containing at least one evanescent operator does not mix 
into physical operators. Together with the statement that evanescent operators do 
not contribute to the matching this proves our method to be consistent at the NLO 
level. As in the case of single insertions one can pick the dki, bki,- ■ ■ completely 
arbitrary and then has to perform a finite renormalization for the double insertions 
containing an evanescent operator in flSSp . This statement remains valid also in 
higher orders of the SU(N) interaction, which can be proven analogously to the 
proof given by Dugan and Grinstein || for the case of single insertions. 

Now we use the findings above to show the nonvanishing terms in (fi7|) explic- 
itly for the physical submatrix: 



J** 

I knl 



4 

-2 



Zi 



1,(2) 



kn,l 
1,(1)1 



kE lk 



Z -1,(D 
7-1.(1) 



kn,E-, 



-1,(1) 



E,un,l 



z\ 



E u ,l 
-1,(1) 



nE ln' 



-1,(1) 



kE ln ,,l 



(49) 
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The last equation encodes the following rule for the correct treatment of evanes- 
cent operators in NLO calculations: The correct contribution of evanescent opera- 
tors to the NLO physical anomalous dimension tensor is obtained by inserting the 
evanescent one-loop counterterms with a factor of | instead of 1 into the coun- 
terterm graphs. Hence the finding of || for a single operator insertion generalizes 
to Green's functions with double insertions. Here the second term in (49) cor- 
responds to the graphs with the insertion of a local evanescent counterterm into 
the graphs depicted in fig. [I], while the last to terms correspond to the diagrams 
of fig. |2| with one physical and one evanescent operator. 



5.3 Double Insertions: Evanescent Scheme Dependences 

In this section we will answer questions |3| and || from page [14]. Let us first look 
at the variation of the anomalous dimension tensor jkk',! on the coefficients a rs . 
First one notices, that the LO 7^, ; is independent of the choice of the a rs . In the 
NLO case one derives in a way completely analogously to the procedure presented 
in section |] the following relation 

7$,, (a') = 7& (a) + [D ■ (a' - a)] ks 7$, + [D • (a' - a)] k , s 7 £ (50) 



with the diagonal matrix D from (p0|) . Note that the indices only run over the 
physical subspace. 

The variation of the anomalous dimension tensor jkk',i with the coefficients 
a rs again vanishes in LO, in NLO we find the transformation 



(i) t~\ 1 (o) 
Tkk'i W + % k \, 



+ 



(°), 

kj 



1,(1) 



(0) 

f 



Zi 



iEi 
1,(1) 



fa' — a] 



2f3o Z[ 



-1,(1) 



kk',E u 



\a — a] 



if, Ei 



\a! - a] 



kk',E u 



~ 1 

a 



~(o) 

■J si 



(51) 



As in the case of single insertions, up to the NLO level there exists no de- 
pendence of 7 on the coefficients b rs and also no one on the b rs . This provides a 
nontrivial check of the treatment of evanescent operators in a practical calcula- 
tion, when the b rs , b rs are kept arbitrary: the individual renormalization factors 
Z each exhibit a dependence on the coefficients b rs , b rs but all this dependence 
cancels, when the Z's get combined to 7. 

Next we will elaborate on the scheme independence of RG improved physical 
observables. First look at the solution of the inhomogeneous RG equation ([43|) 
for the local operator's Wilson coefficient: 



Ci(g^)) 



UP (g(ri,g )Q,(g ) 

2 / N I 9 ^ 



+ 



Oil' + : „ 9 J W 



dg'U^(g^),g') 



90 



,/2 



>kk' 



167T 2 



J> 



kk' 
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,'2 



$nn' r, J nn' 

l07T z 



,'2 



&mm' 



167T 2 



7 J mm! 



U$t (9', 9o) 

uj£) v {g', go) 



,/2 



8i±i — 



1 

lnm,k' L 



Pl (0) _ 1nm,k' 
02^ nm,k' n 

Po Po 



16tt 2 






11 


g 


167T 2 






9 < 1 




167T 2 


\ 



Cv(go 



C v '{g ) 
(52) 



Here the matrices U^°\ denote the LO evolution matrices stemming from 
the solution of the homogeneous RG equations for the Wilson coefficients C, C, 
which reads 



y(0) J 



U {0) (g,go) 



go 
g 



(53) 



We have not labeled the evolution matrices with the renormalization scales //, /xq 
but rather with the corresponding coupling constants g (/i) and g Q = g (/i ). The 
matrix J is a solution of the matrix equation [TC 



J + 



7 



(0) J 



2(3 



J 



7 



(i) 



T 



Wo ' Wl 



+ 



77 



{oy 



(54) 



The matrices J are defined analogously in terms of 7. If 7 transforms 

according to fl3~T|), we know from 



TOU that J transforms as 



J(a') = J(a) - [D ■ (a' - a)f 



(55) 



which can be easily verified from (|32"D. Hence after inserting (^), (0) and (|55|) 
into (^) one finds the independence of C\ from the coefficients a^i- 

In a way similar to the one described above, one treats the scheme dependence 
coming from the coefficients a^. Here some work has been necessary to prove 
the cancellation of the scheme dependence connected to g' 2 Jkk' and ^ mk t in flS"^): 
Although it is not possible to perform the integration in fl52D without transforming 
some of the operators to the diagonal basis, one can do the integral for the scheme 
dependent part of (0), because the part of the integrand depending on a^'s is 
a total derivative with respect to g. There is one important difference compared 
to the case of the dependence on the a^'s: A scheme dependence of the Wilson 
coefficient stemming from the lower end of the RG evolution remains. This is 
a well-known feature of RG improved perturbation theory [|K|. This residual 
dki dependence must be canceled by a corresponding one in the hadronic matrix 
element. If the matrix elements are obtained in perturbation theory, one can 
show that the dependence of the C\ gets completely resolved. Finally, as in 
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the case of single insertions |ID |, one can define a scheme- independent Wilson 
coefficient for the local operator 



cm 



l07T Z 



Cv + ^ • r nmjl ■ Cf (/i) Cg>00 + O 



(56) 



which multiplies a scheme independent matrix element defined accordingly. It 
contains the analogue of f in pi for the double insertion 



2 J QnQr 



(0) 



' nm,l \ I . 



167T 



(57) 



6 Inclusive Decays 

Inclusive decays are calculated either by calculating the renormalized amplitude 
and performing a subsequent phase space integration and a summation over final 
polarizations etc. (referenced as method 1) or by use of the optical theorem, which 
corresponds to taking the imaginary part of the self-energy diagram depicted in 
fig. |3| (method 2). This figure shows that inclusive decays are in fact related to 




Figure 3: The lowest order self-energy diagram needed for the calculation of 
inclusive decays via the optical theorem (method 2). 

double insertions, but in contrast to the case of section [5] they do not involve local 
four-quark operators as counterterms for double insertions. In fact, even local 
two-quark operator counterterms would only be needed to renormalize the real 
part, but the imaginary part of their matrix elements clearly vanishes. The only 
scheme dependence to be discussed is therefore the one associated with the a^'s, 
friz's, etc., as there are no a^s b^s, etc. involved. 

To discuss the dependence on the a^'s it is nevertheless advantageous to con- 
sider method 1, i.e. the calculation of the amplitude plus the subsequent phase 
space integration. From section ^ we already know most of the properties of 
the RG improved amplitude: At the upper renormalization scale the properly 
renormalized evanescent operators do not contribute and the scheme dependence 
cancels. Further we know the scheme dependence of the (RG improved) Wil- 
son coefficients at the lower renormalization scale, because with ( |3"3"D and (pi|) 
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we can use the result of [jTU| . What we are left with is the calculation of the 
properly renormalized operators in perturbation theory, i.e. with on-shell exter- 
nal momenta. Clearly the form of the external states does not affect the scheme 
dependent terms of the matrix elements, they are again given by (|32|) and there- 
fore trivially cancel between the Wilson coefficients and the matrix elements, 
because the scheme dependent terms are independent of the external momenta. 
Since we now have a finite amplitude which is scheme independent, we may con- 
tinue the calculation in four dimensions and therefore forget about the evanescent 
operators. The remaining phase space integration and summation over final po- 
larizations does not introduce any new scheme dependence, therefore we end up 
with a rate independent of the a^'s, few's, etc. 

Alternatively one may use the approach via the optical theorem (method 2). 
Then one has to calculate the imaginary parts of the diagram in fig. |3| plus gluonic 
corrections. Of course the properly renormalized operators have to be plugged 
in: 

Im (6^ n 6[ cn ) (58) 

One immediately ends up with a finite rate. What we only have to show is the 
consistency of the optical theorem with the presence of evanescent operators and 
with their arbitrary definition proposed in (|5|), (|i6|). This means that evanescent 
operators must not contribute to the rate, i.e. diagrams containing an insertion 
of one or two evanescent operators must be of order e 

Im (% [O k y en 6p) = O (e) and Im (£< [O k } Icn E 3 [O^) = O (e) . (59) 

As in the previous sections one can discuss tensor integrals and Dirac algebra 
separately leading to fl59|) . 

7 Conclusions 

In this work we have analyzed the effect of different definitions of evanescent oper- 
ators. We have shown that one may arbitrarily redefine any evanescent operator 
by (D — 4) times any physical operator without affecting the block-triangular 
form of the anomalous dimension matrix, which ensures that properly renormal- 
ized evanescent operators do not mix into physical ones. Especially one is not 
forced to use the definition of the evanescent operators proposed in [|J , whose im- 
plementation is quite cumbersome. Then we have analyzed the renormalization 
scheme dependence associated with the redefinition transformation in the next-to- 
leading order in renormalization group improved perturbation theory. We stress 

3 ^We discard problems due to infrared singularities and the Bloch-Nordsiek theorem. At least 
in NLO one can use a gluon mass, because no three-gluon vertex contributes to the relevant 
diagrams 
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that it is meaningless to give some anomalous dimension matrix or some Wilson 
coefficients beyond leading logarithms without specifying the definition of the 
evanescent operators used during the calculation. In physical observables, how- 
ever, this renormalization scheme dependence cancels between Wilson coefficients 
and the anomalous dimension matrix. One may take advantage of this feature by 
defining the evanescent operators such as to achieve a simple form for the anoma- 
lous dimension matrix. Then we have extended the work of || and || to the 
case of Green's functions with two operator insertions and have also analyzed the 
abovementioned renormalization scheme dependence. For this we have set up the 
NLO renormalization group formalism for four-quark Green's functions with two 
operator insertions, we have derived the renormalization scheme dependence of 
the corresponding anomalous dimension tensors and defined scheme-independent 
Wilson coefficients. Finally we have analyzed inclusive decay rates. 
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